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We establish a relation between the Raman response function in the Big channel and the electronic contri¬ 
bution to the nematic susceptibility within the spin-driven approach to electron nematicity of the iron based 
superconductors. The spin-driven nematic phase, characterized by the broken C 4 symmetry, but unbroken 0(3) 
spin-rotational symmetry, is generated by the presence of magnetic fluctuations associated with the striped 
phase. It occurs as a separate phase between Tm and Ts in systems where the structural and magnetic phase 
transitions are separated. Detecting the presence of nematic degrees of freedom in iron-based superconduc¬ 
tors is a difficult task, since it involves measuring higher order spin correlation functions. We show that the 
nematic degrees of freedom manifest themselves in the experimentally measurable Raman response function. 

We calculate the Raman response function in the tetragonal phase in the large N limit by considering higher- 
order Aslamazov-Larkin type of diagrams. They are characterized by a series of inserted quartic paramagnon 
couplings mediated by electronic excitations that resemble the nematic coupling constant of the theory. These 
diagrams effectively account for collisions between spin fluctuations. By summing an infinite number of such 
higher order diagrams, we demonstrate that the electronic Raman response function shows a clear maximum at 
the structural phase transition in the Big channel. Hence, the Raman response function can be used to probe 
nematic degrees of freedom. 

PACS numbers: 74.25.nd, 74.70.Xa, 74.20.Mn, 74.25.Ha 


I. INTRODUCTION 

Iron-based superconductors show rich phase diagrams, with 
the high-temperature superconducting dome being in the close 
proximity to an antiferromagnetic striped phased that sets in 
at a temperature Tm- In addition, a structural phase transi¬ 
tion at Ts, from the high-temperature tetragonal phase into 
an orthorhombic phase, has been shown to closely follow the 
magnetic transition^r.^, i.e.: Tg > Tm ■ It was proposed 
that spin-fluctuations, associated with the striped phase, lead 
to emergent electronic nematic degrees of freedom at higher 
temperaturesi^i^ These electronic nematic degrees of free¬ 
dom then couple to the lattice and induce the structural phase 
transition to the ortorhombic phase 

There is a mounting evidence for an electronic ne¬ 
matic state: resistivity-anisotropy measurement s * and 
the measurement of the elastoresistancer^ the observed 
anisotropies in thermopower^^ optical conductivityr^i^ 
torque magnetometry^ and in STM measurementsMea¬ 
surements of the elastic constants showed that the shear 
modulus strongly softens in the high temperature tetrago¬ 
nal phaseii2i22ri21 A theoretical analysisi^ based upon nematic 
fluctuations due to a strong magneto-elastic coupling showed 
that the inverse shear modulus is proportional to the suscep¬ 
tibility of the nematic order parameter Xnem, which diverges 
at the structural phase transition, explaining its softening. The 
most direct evidence for the magnetic origin of nematicity so 
far is the scaling of the shear modulus and the NMR spin- 
lattice relaxation rate, seen in iron-pnictides3 An interesting 
open issue in this context is the lack of such scaling behavior 


in iron-chalcagonidesi^ 

A relation between nematicity and the Raman response of 
iron based superconductors was already studied in Ref. 
where the Kramers-Kronig transform of the Raman response 
was compared with the shear modulus, as well as in Ref. 

Here, we demonstrate, based on an explicit micro¬ 
scopic theory that i) there is no enhancement of the electronic 
Raman response function in the i? 2 g channel upon lowering 
the temperature, ii) that the Raman response function devel¬ 
ops a pronounced peak at the structural phase transition in the 
Big symmetry, and iii) that there is some response in the Aig 
channel, which weakens as the temperature is lowered towards 
the structural transition temperature. 

We start from the spin-driven scenario for the nematic 
phase, in which magnetic fluctuations stabilize a nematic 
phase, characterized by the broken C 4 symmetry. The Raman 
response function measures the electronic density-density cor¬ 
relator weighted by appropriate form factors. Since electrons 
interact with spin fluctuations, the latter will manifest them¬ 
selves in the Raman response function in the form of correc¬ 
tions to the electron self energy and the Raman vertex, for¬ 
mally expressed in terms of Aslamazov-Larkin diagrams<2S 
We show that the leading order Aslamazov-Larkin (AL) di¬ 
agram supports only the Aig and the Big symmetry, but not 
the B 2 g symmetry, which explains the lack of enhancement 
of the Raman response signal as one approaches the struc¬ 
tural transition in the B 2 g channel, as seen in experiments 
However, this leading order approach cannot account for the 
rapid increase in the amplitude of the Raman response func¬ 
tion in the Big channel, as one approaches the structural tran¬ 
sition, as seen in the experiments of Refs. mui Instead it 
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would predict a similar increase only at the magnetic phase 
transition. Therefore, we go beyond the leading order ap¬ 
proximation, and take into account collisions between spin 
fluctuations that become more and more important as one ap¬ 
proaches the nematic / structural transition. Our approach is 
based on the exact same collisions between spin-fluctuations 
that led to the emergence of spin-induced nematicity in the 
first place. Formally this is accomplished by inserting a se¬ 
ries of quartic paramagnon couplings, mediated by electronic 
excitations, into the Raman response function. Such quartic 
couplings contain a product of four fermionic Green’s func¬ 
tions and give rise to a peak of the electronic Raman response 
function at the structural phase transition in the Big channel. 
On the other hand, if we re-sum such higher order AL dia¬ 
grams in the Aig channel, this will lead to the suppression of 
the Raman response in the aforementioned channel. 

Here we demonstrate that the low frequency Raman re¬ 
sponse in the Big channel is given by 


Rbi, (w) 


Ro{^) 


( 1 ) 


where w is the frequency difference between incoming and 
outgoing photons and g the multi-index for momentum and 
frequency. stands for the leading order Aslamazov- 

Larkin diagram, Xg is the magnetic susceptibility, and g the 
nematic coupling constant of the theory. On the other hand, 
the susceptibility of the nematic order parameter of our model, 
in the large N limit is given by 


Xnem — 


IqX 


2 

Q 


1 5stat Jg Xq 


( 2 ) 


where in a purely electronic theory g = pstat- In a purely 
electronic theory, this would then lead to the divergence of 
the Raman response function at the structural phase transi¬ 
tion. However, one needs to include the effect of the lattice 
dynamics^ in order to analyze this problem. We do so by in¬ 
troducing nemato-elastic coupling and And that, in this case, 
2 

ffstat = 5 + ^ is shifted^ Here 7 ei is the elasto-nematic 

coupling constant, and Cg the bare value of the orthorhombic 
elastic constant. We show that when magnetic and structural 
phase transitions are split^ii^ this leads to a maximum of the 
amplitude of the electronic Raman response function in the 
Big channel at the structural phase transition, in agreement 
with the recent experiments^ The Raman response function 
could then be used to probe the dynamic excitation spectrum 
of the nematic degrees of freedom, similar to inelastic neutron 
scattering that probes the dynamic spin excitation spectrum. 

The paper is organized as follows. In Sec. |II]we present the 
microscopic model for the spin-driven nematic phase. We cal¬ 
culate the effective action and analyze it in the large-A^ limit, 
where N'^ — 1 is the number of components of the collective 
paramagnon held. Following Ref. 0, we derive the condi¬ 
tion for the susceptibility of the nematic order parameter to 
diverge. In Sec. |III]we then show how to calculate the Ra¬ 
man response function using a diagrammatic approach. We 
first calculate the leading order Aslamazov Larkin diagram. 


and show that there is no response in the i? 2 g channel, and 
a finite response in the Big and the Aig channels. We then 
calculate higher order diagrams that take into account colli¬ 
sions between spin-fluctuations. Finally, after summing an 
infinite number of these higher-order diagrams within a con¬ 
trolled 1 /N expansion, we show i) that the maximum of the 
Raman response function in the Big channel occurs when the 
nematic susceptibility diverges, i.e. at the structural phase 
transition, and ii) that the amplitude of the Raman response 
function in the Aig response gets suppressed. We present our 
conclusions in Sec. EYl 


II. MICROSCOPIC MODEL: SPIN DRIVEN NEMATICITY 


Two different approaches have been proposed in order to 
explain the origin of nematic phase in pnictides and its re¬ 
lation to the magnetic phase - the orbital scenario^— and 
the spin-driven nematic scenarioi^!^ For a discussion of these 
approaches see for example Ref. Here we follow the ap¬ 
proach of a spin-driven nematic state. In this scenario, the 
nematic phase is stabilized by magnetic fluctuations that are 
associated with the stripe density wave (SOW) phase. The or¬ 
der parameter of the SDW stated can be characterized by an 
0(3) X Z 2 manifold^^ - 0(3) is the spin-rotational sym¬ 
metry and Z 2 a discrete symmetry associated with the choice 
of the ordering wave-vector, Qx = (tt, 0) or Qy = (0,7r). 
Let the two order parameters associated with these two or¬ 
dering wave vectors be Ax and Ay respectively. The SDW 
state is characterized by broken 0(3) and Z 2 symmetries. On 
the mean-field level the breaking of Z 2 and 0(3) symmetry 
occurs simultaneously. However, when one includes fluctua¬ 
tions, these transitions can be split. In case of joint transitions, 
they are usually both first order transitions The criterion for 
breaking the discrete Z 2 symmetry via a second order tran¬ 
sition is a threshold value of the magnetic correlation length 

Decreasing the temperature leads to an increase of Be¬ 
fore the correlation length diverges at the magnetic phase tran¬ 
sition temperature, the threshold value will be reached and 
spin-driven nematicity sets in. This naturally explains why 
the magnetic and structural phase boundaries are correlated 
and leads to an intermediate phase with Z 2 symmetry break¬ 
ing without 0(3) symmetry breaking. This intermediate state 
is the nematic phase in the pnictides. It is characterized by 
unequal strength of the magnetic fluctuations associated with 
the ordering wave vectors Qx and Qx : ( A^ - A?.) ^ 0, 
but no long range magnetic order, ( Ax,f) = 0 . 

In what follows we will summarize the steps of Ref. 0 and 
outline the mathematical model the for spin-driven nematic 
phase. We start from a simplified itinerant model where we in¬ 
clude the bands near the F—point and the X— and Y— points 
in the Brillouin zone. For our main result no explicit knowl¬ 
edge of the detailed parametrization of the band structure is 
necessary, except for the fact that the band-structure is not 
perfectly nested. However, in order to obtain emlicit numeri¬ 
cal results we use the simplified model of Ref. □. We consider 
parabolic dispersions with 
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FIG. 1; Band structure: the model consists of the central hole-like F 
band, and the electron-like X and Y bands, shifted by Qx = (7r,0) 
and Qy = (0, tt), respectively. 
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where mi are the band masses, cq is the offset energy, and 
fi denotes the chemical potential. The corresponding Fermi 
surfaces are shown in Fig[T] 

In order to study the established stripe magnetic phase, we 
consider the Hamiltonian that contains the interactions in the 
spin channel with momenta near Qx and Qy: 


% = 'Hq + Hint, 

Ho = ^ ' Cz,kCi,kQ,Ci,kQ;, 


2,k 


Hint — ^ ^ 


-q • 


(4) 


Here, c| is the creation operator of an electron with mo¬ 
mentum k, spin a and in the band i. The spin operator is 
given by 


~ '^r,k-f.qct'^“/3'^bk^, (5) 

k 

where \ap denotes the iV^ — 1 component vector of the gen¬ 
erators of the SU{N) algebra. In the case TV = 2 it holds 
•^ 0/3 = f fa/S with vector of the Pauli matrixes cr. Ug is the 
coupling in the spin channel, which can be expressed in terms 
of density-density and pair-hopping interactions between hole 
and electron pockets 


The partition function is given by 

Z = y (6) 

where [3 = T~^ is the inverse temperature. Since, the in¬ 
teraction Hamiltonian is quadratic in the fermionic spin q, 
we can decouple it using a Hubbard-Stratonovich decoupling 
in the spin-channel. This way, we eliminate the quartic in¬ 
teraction between fermions at the expense of a functional in¬ 
tegral over two additional bosonic fields Ax and Ay, with 
N"^ — 1 components. The bosonic fields couple linearly to the 
corresponding fermionic spin densities. After introducing the 
spinor 


'^k ~ ( Cr,ka '^X,ka '^y.ka ) ’ 

where a denotes every possible value of the SU{N) spin in¬ 
dex, we can write the partition function as: 

Z = [ (8) 


with the action: 

^ [vk, A,] = - / + ^ f + ^y) ■ (9) 

Jk Jx 


Here, the matrix of the inverse Green’s function Q/. ^ is given 
by: 

Ga^ = Go,1-Va, ( 10 ) 

with the bare term: 

/Gr.fc 0 0 \ 

ao,fe= 0 Gx.fe 0 , (11) 

\ 0 0 Gy,fc / 


and the interacting term: 



( 12 ) 


Gi,fe = Gi,k3 with G^l = iuJn — Ci.k and N x N unit matrix 
1. We invert the matrix equation (fTOl) by expanding the geo¬ 
metric series and obtain the following expression for Qa that 
we will use later-on: 


Sa = ^(eoVA)"eo. (13) 

n—0 


A. Effective action in the large-A expansion 

In this section, we first show how to obtain the Ginzburg- 
Landau expansion of the effective action in powers of the spin 
fluctuation fields Ax,y in the limit of large N, in the spirit 
similar to that of Ref. □ where only N = 2 was considered. 
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Next, we re-formulate this effective action in terms of the col¬ 
lective nematic Ising variable (p, and analyze the equation of 
state for p. We deduce the condition for the onset of the ne¬ 
matic phase by examining the susceptibility of the nematic 
order parameter. We begin by integrating out the fermionic 
degrees of freedom from Eq. It follows: 

Z = (14) 

with action: 

5'eff[Ax,A'^] = —Trln (1 — C/oVa) 

+ -/(Ax + A^). (15) 

UsJx 

Here, Tr (• • •) refers to sum over momentum, frequency, spin, 
and band indices. We further expand in powers of Ax,y to 
obtain: 

5eff[Ax,Ay] = iTr(^o.fcVA)" + iTr(eo.feVA)^ 

+ — [ (Aj^ + Alr) . (16) 

J X 


with coefficients 

2 1 f 

’"0 =-^ o / ^x,kGr,k, 

Us 2 J 

u = J Gr,fc(G'x,fe + Gr.fc)^, 

9 = J^GlkiGx,k-GY,k)^. (21) 

To, u and g have been calculated as a function of temperature 
and band parameters in Ref. 0. It was found that u > 0 
and u > 5 in general. The coupling g vanishes for circular 
electron pockets, but is positive for a non-zero ellipticity. 


B. Nematic susceptibility in the large-N expansion 

In order to investigate the possibility of the nematic tran¬ 
sition occuring before the magnetic transition, we follow 
the steps of Ref. 0, and introduce two auxiliarly Hubbard- 
Stratonovich scalar fields (p and ip to decouple the quartic 
terms in the action (l 20 l i: cp —>■ A^^ — Ay and ip Aj^ + Ay- 
The resulting effective action is given by 


After using a series of identities for the generators of the 
SU{N) algebra, needed to evaluate the above traces (for de¬ 
tails see lA 21 i. we arrive at the following effective action in the 
large N limit: 

SeS [Ax, Ay] = ^ ro^iAf + UijAfA^, (17) 

1 i,j 

with the coefficients: 

2 If 

1^0,1 = -1" n / Gr,kGi k, 

Us 2 Jk 

Uij = J^Gr,kGi,kGj^k- (18) 

We used the notation = Tj^f 7^^- The index k = 

(k, w„) combines the momentum k and the Matsubara fre¬ 
quency ujn = {2n -f 1) ttT. 

Using the identities 


SeS 



Jx\2g 2uj 


f V' (Ax + Ay) + f {(p + hn) (Ax - Ay) , 
J X J X 


( 22 ) 


and we have added a field hn conjugate to the nematic or¬ 
der parameter — Ay. This term is needed in order to 
calculate the susceptibility of the nematic order parameter. 
A finite value of p implies non-zero expectation value of 
^ = (A^ — Ay) ^ 0 and the system develops nematic 
order. The large-A mean field value of ip is always non-zero 
and describes the strength of magnetic fluctuations. In case 
of split magnetic and structural phase transitions, there is no 
magnetic order right below the structural transition temper¬ 
ature, i.e (Ax,y) = 0. Next we integrate out the — 1 
component fields Ax.y- If we further rescale the coupling 
constants to g = g{N‘^ — 1 ) and u = u{N‘^ — 1 ), required 
to reach a sensible large-A limit, the effective action can be 
written as 


[ Gr.kGx,k = [ GT,kGY,k, 
Jk Jk 


/^2 /^2 _ 


/^2 r^2 


(19) 


valid because the underlying Hamiltonian obeys the full C 4 
symmetry, we can write the action in the more convenient 
form: 


>S'eff[Ax, Ay] — ro(Ax + Ay)-f —(A 

-f(Ai-A^)^ 


A|.)2 

( 20 ) 


Ses[lp,p] 



Ip'^ 

2u 






(23) 


We note that the effective action (l23t has an overall pre-factor 
A^. For A ^ 1 the integral over the fields p and p can be 
performed via the saddle-point method, i.e. by analyzing the 
extremum of the action. After solving for dSeS [p, p] /dp = 
dSeS [p, p] /dp = 0 , we obtain the equations of state for p 
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and iIj: 

^ _ Xq^ +'^ _ 

“ {Xq^ +'(P)^ - (4'+ hnf 

i = [ _ t±]hL _. (24) 

^ ‘I {Xq ^ + V”) “('/’ + ^n) 

By differentiating the second equation in (l24l i with respect to 
the conjugate field, we find that, for small (p 


94 ^ 9 Ik xl 

dhn 1-ff/fcXr 


(25) 


where, from now on, we have shifted Xk^ ^ Xfc ^ which 
simply corresponds to the re-normalisation of the mass term 
due to fluctuations. Similarly to the result of Ref. liH we 
find that the electronic contribution to the susceptibility of the 
nematic order parameter — Ay is given by 


X 


el 

nem 


/fcXfc 

1 -ff/fcxr 


where Xq ^ is the inverse magnetic susceptibility, and 


(26) 


N f 
^ " 16 Jk 

is the nematic coupling constant of the theory. In Ref. 0 it 
was found that for the classical phase transition in d = 2 
and u/g > 2 the nematic transition pre-empts the magnetic 
transition, i.e. the transition lines are split. Also, the ne¬ 
matic transition was found to be of second order. This is the 
regime we are interested in. What we have calculated so far is 
the purely electronic contribution to the nematic susceptibil¬ 
ity. One, however needs to include the effect of the lattice, as 
was pointed out in Ref. l2lll^ The nemato-elastic coupling is 
given by the following Hamiltonian 


where 


dstat — 9 + 


7el 


(31) 


III. RAMAN RESPONSE FUNCTION 


Raman scattering is a valuable tool to study strongly cor¬ 
related electronic systems^, since it probes lattice, spin and 
electronic degrees of freedom. It has been used to extract in¬ 
formation about the momentum structure and symmetry of the 
excitations in the cuprates^M^.^ and pnictides. The differen¬ 
tial photon scattering cross section in Raman spectroscopy is 
directly proportional to the structure factor S: 

Sq = -[1-b n(a;)] Imiig, (32) 

TT 

which is related to the imaginary part to the Raman response 
function R through the fluctuation-dissipation theorem4^ 
Here, n(w) is the Bose-Einstein distribution function, and 
q = (q, w). Since the momentum oflight is much smaller than 
the typical lattice momentum, one normally replaces q Ri 0 in 
Eq. 

The Raman response functions measures correlations be¬ 
tween “effective charge density” fluctuations p, 

/.l/T 

R{u}) = / dre““^(p(r)p(0)). (33) 

Jo 

The effective density, weighted by the form factors that can 
be changed via the geometry of the photon polarization, is 
defined as 


Pk= Y. Xk>clk+k>,aCi,k',a- (34) 


^nem — Tel 


dr(f){r) {dxUx - dyUy ), 


(28) 


where xa is the nemato-elastic coupling constant and u = 
{ux,Uy) the phonon displacement held. The phonons renor¬ 
malize the nematic coupling constant to a frequency and mo¬ 
mentum dependent coupling 

„2 

g{q,u}) =g + (29) 

where c® is the elastic constant and q the momentum along 
the soft directions. In particular, if one wants to determine the 
location of the nematic phase transition, which is dictated by 
the condition of divergent nematic susceptibility, one needs to 
look at the static limit of the coupling constant, i.e. the limit 

where w is set to zero. This leads to p^tat = 5 -f 4- The full 
nematic susceptibility, including the effect of the coupling to 
the lattice, in the large N expansion is therefore given by 


IkXl 

1 - 5stat Ik xl ’ 


( 30 ) 


a is the spin index, i the band index, and the operator cT ^ 
creates an electron with spin a and momentum k in band i, 
where i = X,Y, E. The function 7 k is related to the inci¬ 
dent and scattered photon polarization vectors and depends 
on the curvature of the bands4^ Here, we will consider the 
intra-orbital contributions to Xk only, as this is the dominant 
process. The multi-orbital nature of different bands has been 
pointed out in Ref. 

In order to determine the Raman response function, we cou¬ 
ple an external source held to the weighted densities and in¬ 
troduce the generating functional Wh according to; 




(35) 


where S [fll, A^] is given in Eq. (|9l). The elements of the ma¬ 
trix Vh in momentum/frequency, spin and band space are 

'^h,k\k2(7a'ij — — A:27ki ; 


(36) 
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with h being the field conjugate to the effective density. The 
Raman response function (l33T l is obtained by differentiating 
the generating functional Wh (l35]) with respect to the conju¬ 
gate field h: 


Rq 


S^Wh 

ShqSh-q 


h=0 


(37) 



Due to the single particle character of the source term, the 
generating functional Eq. dTSl) can be written in the form 

Wh = ^ J 

= Go"-VA-Vh. (38) 

Since Wh contains the action that is quadratic in fermions, we 
integrate out the fermions and obtain: 

Wh = ^ j 

Sh[A,] = 1 / (A^ + A^)-Trlnfc)- (39) 

J X ^ ' 

We further expand: 


Tr In I 


=Trln( 0 -E 

n—1 

Then, using (l3^ and (ITTI i. 


Tr {gAVhY 


(40) 


ShqSh-q 


exp 


-TriGAVh)-^Tr{GAVhf 


h—0 

(41) 


Here S'eff[Ax,Ay] = is the effective action 

given by (fTbl) . We define the matrix 


r« = 


Shq 


(42) 


FIG. 2; Left: Contribution to the Raman response function that con¬ 
tains the self-energy correction to the fermionic propagator. Right: 
A diagram that contains a vertex renormalization correction. 


In order to get the vertex correction, we replace both Qa in 
( l43b by Qa Go'^aGo, which comes from the perturbative 
expansion of Eq. O: 


rV = - 


(^oVa^oT)^ 


• (45) 


In order to get the self-energy correction, we replace one 
of Ga in ( |4T1 i by Ga Go, the other one by Ga 
{GoVaY Go^o get 


= 






SoT [GoVaY ^oT 


(46) 


Due to the integral over the square of the 7 ^ factor, the self¬ 
energy and vertex corrections occur in all symmetry channels. 
If one evaluates the sum + R^ explicitly, in the hot-spot 
approximation, one finds that there are partial cancellation in 
the Aig and in the Big channels, and no cancellations in the 
B 2 g channel. One can easily show that in d = 2 

R^ + R^ (X f -— g oc log^, (47) 

7q ro + 

where we have used tq = where ^ is the magnetic corre¬ 
lation length. 


B. Leading order Aslamazov-Larkin diagrams 


A. Self-energy and vertex correction diagrams 


Next, we analyse the leading order contributions to the Ra¬ 
man response function. These arise from the self-energy and 
vertex correction diagrams depicted in Eig. |2] Both of these 
diagrams arise from differentiating the second term in the ex¬ 
ponential (I 4 TI 1 twice with respect to h 


R 


v,s 

q 


1 


dA,e"^-fdAx.Av]Tj. 


(^Ar)' 


(43) 


The Aslamazov-Larkin contribution to the Raman response 
function, analyzed in Ref. IH arises from differentiating the 
first term inside the exponential in dril l twice, and from re¬ 
placing Ga {Go^aY Go, which comes from the perturba¬ 
tive expansion of Eq. dS): 


Rg = - 


Tr((aoVA)"eor) 


(48) 


and we replace SeS So, where So is the quadratic action 
given by 

So [Ai] = — f (Ax + Ay) 4- —Tr {Go'^a) ■ (44) 


Here S'eff[Ax,Ay] = S'/i[Ai]|?i=o is the effective action 
given by dl^ . 

As we will see below, the key assumption of a description 
based on the Aslamasov-Larkin diagrams is that one neglects 
the interactions between spin fluctuations. In other words, one 
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by Eq. ( |49] | and (fSOl l. only supports the Big and the Aig sym¬ 
metry channels. Let us consider the structure of the terms in 
(|49] | which arise from 


i?, 


( 11 ) 


(c^):= tY. 


2 =X,y,r 


xx(q,f^«)x(q,^n - w). (51) 

The term ( ISTT i can be re-written in the following form 


Rq^\^)= / f /7k7px(q,f^«)x(q,f^n-w) 

^ n -Iq-lkJp 

X (w, ^n, k, q)E)Aig (w, p, q) 
+-E'Big k, q)E)Big (w, p, q)), 

(52) 


FIG. 3: Leading order Aslamazov-Larkin diagram. Raman vertices 
are denoted by black squares. Imaginary part of the Raman response 
function as a function of frequency Im7?o(u;), in d = 2. 


approximates the effective action in ( l48l l by quadratic action 
Eq. (l44l i. While this assumption is frequently justified, it is 
not allowed in the theory of spin-driven nematicity, as we will 
show later. 

The leading order Aslamazov-Larkin diagram, depicted in 
Eig. [3can be calculated as 


Ro{u}) 


i=X,Y,n'^^ 


(49) 


with 


Ai(q, 

Af^(q,f2,w) 


Ap^(q,Ll,w) -F Af^(-q,-fl,-w), 

tY f 7 kGr(k, r-n - w)Gr(k, !/„) 

^ k 


xGi(k-q, -n), 

tY f 7 kGi(k, z/„ - w)Gi(k, i/„) 

n “'k 

xGr(k-q, z/„-n), (50) 


similar to what was found in Ref. HH 


]. Raman response in different symmetry channels 

In the concept of the pairing symmetry in high-temperature 
superconductors successful theoretical models supported by 
experiments have been developed in order to explain the sym¬ 
metry sensitivity of the Raman response function42 Simi¬ 
larly, here, before we explicitly evaluate the leading order 
Aslamazov-Larkin diagram, we analyze the contribution to it 
in the various symmetry channels. Higher order corrections 
that will be discussed later do not alter this symmetry-based 
analysis. We show that the Aslamazov-Larkin diagram, given 


where we have classified the appropriate combinations of 
Green’s functions according to their symmetry into 

q) = T^Gr(k, -w)Gr(k, z/„) 

m 

xG(+^(k - q, Pm - H„), 
q) = T^Gr(k, -w)Gr(k, z/„) 

m 

XG^ ^ (k q, Vjji , 

(53) 

and we have defined G*-^^ = Gx ± Gy- Lrom Eq. (l5Tt . we 
see that the response will be non-zero only for 7 factors in the 
Aig or the Big symmetry. Similarly, using the same line of 
arguments, one can show that all other terms in ( |49] | support 
the Aig or the Big symmetry only. We have thus ruled out the 
response in the i? 2 g channel. 


2. Explicit calculation of the leading order Aslamazov-Larkin 
diagram 


The leading order Aslamazov-Larkin diagram has been 
evaluated in Ref. |28l assuming that the main contribution 
comes from the hot-spot regions and that the momenta of the 
fluctuations are peaked around q Ri Qx,y. After the analytic 
continuation to the real frequencies, we found that the imagi¬ 
nary part of the Raman response function, which is a quantity 
of experimental interest, is given by 

de 

Imi?o(w-F *0“'') = / —[n(e) — n(e-F w)] 

J—00 ^ 

X f Im [x^(e,q)] Im [x^(e-Fa;,q)] , 

(54) 


with the spin propagator in the tetragonal phase given by: 


x^(q,^^) 


1 

^■o + q^ - ’ 


(55) 
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where ro tunes the distance from the magnetic transition, see 
Eq. (ITTI i. In d = 2 the q integral in Eq. ( l54b in can be 
performed exactly, which leads to the following expression 

Im [i?o(w+ i 0+)]^^2 = [ de[n(e+)-n(e_)] 

x[i^(e+)-F(e_)], (56) 

with 

Fix) = — farctan— — ^sgn{x)\ . (57) 

X \ X 2 J 

We defined e± = t Flo 12. The plot of the function (l56l l is 
shown in Eig. [3 In particular one can deduce that, in the 
regime where temperature T is the biggest scale, T > ro, 
Ra{oj)d =2 — ^ for small frequencies lo, while the amplitude 

of the Raman response function scales as i?o“(a;)d =2 — ^ 
in this regime. 

In summary, we have shown that the leading order 
Aslamazov-Larkin diagram gives a non-zero response in the 
Big and Aig symmetries only. It predicts the divergence of 
the Raman response at the magnetic transition, and does not 
carry any signatures of the structural transition. We therefore 
need to go beyond the leading order Aslamazov-Larkin dia¬ 
gram. 


C. Higher order Aslamazov-Larkin-like diagrams 


Next, we go beyond the quadratic action approximation for 
SeS in (l48l l. and include the full quartic action to evaluate the 
Raman response function. As we will show, diagrammati- 
cally this corresponds to inserting a series of fermionic boxes 
that resemble the structure of the nematic coupling constant g 
into the leading order Aslamazov-Larkin diagram in the Big 
symmetry. These diagrams take into account the collisions 
between spin fluctuations which were not accounted for in the 
leading order Aslamazov-Larkin diagram. 

Eirst we show how these terms arise from the diagrammatic 
expansion. We start from Eq. (|48]) . but this time we go be¬ 
yond the quadratic approximation for the effective action, and 
include quartic terms 




Tr((goVA)"0or) 


(58) 


where 


5eff[A,] = 5o[A,] + iTr(goVA)", (59) 

with the bare action 

So [Ai] = — f (A^ + Ay) + —Tr(^oVA) • (60) 

tts Jx J 


We further expand the exponential 


g-iTr(GoVA)^ 


E 

m—O 


1 

ml 


^Tr(c;oVA)" 


( 61 ) 



FIG. 4: Re-summed Raman response function. The resummed box 
Bajs is shaded grey. The first index of the matrix B denotes the 
type a = X,Y of entering spin fluctuations, and the second index 
the type of exiting spin fluctuations. We insert the grey shaded box 
into the Raman reponse, and make some further simplifications to 
evaluate the Raman response function (see the main text). 


to obtain 




m—0 


(62) 


where we averaged the following terms with respect to the 
Gaussian collective spin action: 


r(™) — 

J-Lg 


^Tr(c;oVA)" 


Tr((goVA)"0or)]" 


In order to evaluate the expectation values one performs con¬ 
tractions of the A fields. We obtain a series of diagrams that 
look like the leading order Aslamazov-Larkin diagram with 
an arbitrary number of inserted fermionic boxes, depicted in 

Fig- El 

The higher order diagrams effectively take collisions be¬ 
tween spin fluctuations into account, which have been ne¬ 
glected in the leading order Aslamazov-Larkin diagram. As 
one approaches the transition line, collisions between spin 
fluctuations become more and more important and one would 
anticipate significant changes in the Raman response function 
due to these processes. As we will show, the re-summation of 
boxed Aslamazov-Larkin diagrams will lead to the maximum 
of the Raman response function at the structural phase transi¬ 
tion in the Big channel, and the suppression of the response 
in the Aig channel. 

The next task is to re-sum an infinite number of such di¬ 
agrams. Every box can be characterized by two indices: 
the first one denotes the type of incoming spin fluctuations, 
this can be either X or Y and the second one denotes the 
type of exiting spin fluctuation. Let us denote this box Ba/ 3 . 
Summing all boxed diagrams can be most efficiently ex¬ 
pressed as: 


R{lo) — Rq{uj) + / AQ,(a;,n, q) 

X - w) 

X Bap{ci,q',n,n',Lo)x{q.',^') 

X x(q',f^'- w)A;3(w,n',q'). 


(63) 
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For our analysis it is sufficient to calculate the box at mo¬ 
menta q, q' Ki Qx.y and zero frequencies, which is justified 
for small incoming Raman frequency w, and in the classical 
regime relevant near a finite temperature phase transition. We 
write the Raman reponse function in the tetragonal phase: 


R{u}) « i?o(w)-f / Aa(w,0,q)Ba;3 

j q,q' 

xx^(q,0)x^(q',0)A;3(w,0,q'), (64) 

where Rq{ui) is the leading order diagram. 

The symmetry of the fermionic triangle is such that 

_ _\Blg 

; 

A^'® = (65) 

This relation can be obtained by simply performing a coordi¬ 
nate system rotation by 7r/2 inside the momenta integrals in 
(ISOl l. This allows us to explicitly perform the matrix multipli¬ 
cation, which yields: 

^Bigi^^) = Ro{i^) + Roi^^){Bxx — Bxy) ( x^(q, 0 ), 
^Aig(w) = Ro{uj) + Roiu}){Bxx + Bxy) j X^(q, 0). 

Jq 

( 66 ) 

Next we need to determine an expression for the full box 
Bap, i e. perform a sum over the leading box-diagrams within 
the 1/7V expansion. This is illustrated in Fig. |4]and can be 
written as: 


BaP = Bap+BasBsp / X^(q^0) + ■•■ 

Jq' 

oo / c \ — 

= Y.^B^)°‘p[ , (67) 

m=l \7q / 

The matrix B was deduced from Eq. (l63t and Eq. (IB4b . Eor 
details about explicit evaluation of the SU {N) trace pre-factor 
(which arises from contractions of products of A matrices in 
(fT2l i) for boxed diagram containing arbitrary number of boxes 
m, please see Appendix iBl The matrix B of irreducible boxes 
is then given as 


^ ^_N. (dxx gXY\ 
where we used the abbreviation 


gxx 

gxY 


/^2 /^2 


G 2 


( 68 ) 


(69) 


and used that by symmetry: 

The mth power of the matrix B is given by 


B^ 


2 V 8 ; U3+ -5-) 


(sT - 9-) 

{g+ + 9-) 


,(70) 


where g± = gxx ± gxY- Erom this analysis follows that 




= Ro{io) + Ro{uj) 


g/qX^(q,o) 

1 -g/aX^(q,o)' 


Ra,^{u:) = Roiu;) f2 ( / X^(q,0) 

m=0 V ° / V7q 


w/qX^(q,o) 


where 


= Roiu:) + Roiiv) (71) 

1 + ujq xnq,o) 

N f 

^ " 16 A 


is precisely the nematic coupling constant of Eq. (l2ni for the 
effective action, and u is the other quartic term in Eq. (l2Tl i. 
with M > 0, as found in Ref 0 Erom (fTTT i. we see that the 
Raman response in the Aig channel gets suppressed, due to 
the term in the denominator, which grows as one approaches 
the transition. On the other hand, in the Big channel, after 
performing the analytic continuation to real frequencies and 
taking the imaginary part, we get that 


huRB^^iuj) = Im[i?o(u;)] (1 + pXnem) , (73) 


where 

I fqXHq,0) 

Xnem -i ~ r 9/^ n\ 

1 - gjq xnq,o) 


(74) 


is the electronic contribution to the nematic susceptibility cal¬ 
culated in the large N limifl^ for the model described in Sec. 
Ill Al As was pointed out in Ref. [3^ the enhancement of the 
static nematic coupling constant OTI) does not enter the Ra¬ 
man response, due to the fact that the Raman response oper¬ 
ates in the dynamical limit (q = 0 and finite oj), and the static 
and dynamic limits do not commute^. At the nematic / struc¬ 
tural phase transition the nematic susceptibility (l30t diverges, 
and 



(75) 


Consequently the Raman response function in the Big chan¬ 
nel, given by (l73l l. has a maximum rather than a divergence at 
the structural phase transition. 


IV. CONCLUSION 

In summary, we have shown that the Raman scattering can 
be used as a tool to probe the nematic phase in pnictides. We 
have presented a calculation that demonstrates that, in the low- 
frequency limit, and large N limit, the Raman response func¬ 
tion shows a clear maximum at the structural transition tem¬ 
perature in the Big channel. 
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In our model, the electronic nematic phase in pnictides 
is stabilized by spin-fluctuations associated with the striped 
phase, and occurs as a thin sliver above the magnetic transi¬ 
tion temperature. In order to calculate the Raman response 
function, we have gone beyond the leading order Aslamazov- 
Larkin diagram, and included higher order diagrams that con¬ 
tain a series of quartic paramagnon couplings, mediated by 
electronic excitations. Such quartic couplings contain a prod¬ 
uct of four fermionic Green’s functions and include the ef¬ 
fect of collisions between spin fluctuations. When re-summed 
these diagrams lead to the maximum of the electronic Raman 
response function at the structural transition in the Big chan¬ 
nel, and the suppression of the response in the Aig channel. 

The method that we developed analysed the Raman re¬ 
sponse function only in the regime of small frequencies. It 
would be desirable to extend it to the entire frequency range, 
such that one can analyse the entire shape of the Raman re¬ 
sponse function as a function of temperature, and possibly be 
able to extract some information about the dynamical nematic 
susceptibility. 

Further, one might expect a charge driven nematic phase 
to have similar signatures in the Raman response function. 
This could be relevant to the peculiar case of FeSe, where the 
nematic phase has been detected, but no magnetic phase has 
been seeni^iS In order to do so, we would need to develop a 
theoretical method that goes beyond the large N expansion. 
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summed over. 


{Aj, Afc} 

= + djfc/A; djki = dkji, 

(Al) 

[Aj, A/e] 

= ^jkl = —^kjl, 

(A2) 

Aj A/e 

II 

§ "■ 

+ 

tol h- 

(A3) 


:= djki+ifjki- 

(A4) 


Here duji is symmetric under the exchange of its indices, 
while ikji is antisymmetric under the exchange of neighbour¬ 
ing indices. Eurther, some useful relations for the summations 
of structure constants can be derived,i^^ which read 



- 4 . 


^akl^bkl 

II 

(A5) 

^akl ^bkl 

= N6ab, 

(A6) 

O 

II 


(A7) 


I 

Useful identities that involve the traces of the SU{N) ma¬ 
trices are 


Tr(A,) = 0, (A8) 

Tr(l) = W, (A9) 

Tr(A,A,) = (AlO) 

In order to analyse the trace of the product of four SU{N) 
generators we evaluate 


Tr {\i\j\k\i) 


Tr 








(All) 


where we used the identity (1A31) in the first line, as well as 
Eq. (IA8I) and Eq. (lAlOl) in the second line. These results will 
be of importance for the subsequent analysis of higher order 
diagrams. 


2. Effective action from tr log expansion 


Eirst we calculate the quadratic terms in the free energy ex- 
Appendix A: Effective action of the SU(N) fermionic model pansion. This is given by 


1. Some useful SU(N) identities 

Here, we present some useful identities for the structure 
constants of SU{N). They have been used to determine the 
scaling of the boxed Aslamazov-Larkin diagrams with N, and 
to develop the Ginzburg-Landau expansion of the effective ac¬ 
tion in powers of spin fluctuation fields A (see Sec. |II]i. We 
begin by listing some standard SU{N) identities for the ma¬ 
trices Xi, where i = 1, ..N'^ — 1. All repeated indices are 


N^-l 


iTr(goVA)" = E / Gc.fcGr.fe ^ Tr (A,A,) 


bi=i 


= ^E / Gc<,fcGr,fe|A 


2 

a I I 


(A12) 


where a = X^Y and we used the identity (lAlOb . 

Next we calculate the quartic term in the free energy expan¬ 
sion 









11 


^Tr(eoVA)^ = iTr(A,A,AfcAO 

a^^,y 

X ^ SaaA^A^sA^A^, 
a=X.Y 


We note that in the large N approximation there are no A^r ■ 
Ay terms in the action; however if one considers corrections 
to large iV these terms might appear in the effective action. 

Appendix B: Identities containing products of traces of SU{N) 
generators 

In this appendix we derive further identities for the traces 
of the SU(N) generators, which have been used to deduce the 
dependence of the Aslamazov-Larkin boxed diagrams on N. 
In particular, we would like to calculate 


with 


(A13) 


gxx 


gvY 


'^X.k 




Tm := Tr(Ai^AijTr(Ai2AjjAj3AjJ 
X Tr (A 24 A 23 A^g A^q )... 

XTr {^i2m^i2m-l^i2m + l^i2m. + l) 

xTr . (Bl) 


gxY 


gYX = J Gjc.fcGF.fcGpj., 


(A14) We begin by considering to = 1. Written out explicitly, it 
follows: 


and we used the notation a for ’not a’, i.e. if a = AT then 
a = Y and vice versa. We further substitute the identity (lAl lb 
in (IA13b . to write 

^Tr {GVAf = Ki + K 2 , (A15) 

where 

^^2 = E ^R*,pRfc;pALA^„A^AL 

+ ^ (A16) 

a=X,Y 

Since K 2 ^ N~^, while Ki ^ N~^, the term K 2 can be 
omitted in the large N limit. 

Combining ( IA16b and (IA12b . the effective action in the 
large N limit can be written as 

5'eff [Ax, Ay] = ^ ?'o,iA^ + ^ UyA^A^, (A17) 
i i,j 


Ti = Tr(A,Aj)Tr(AfeA/)Tr(A,A,AiAfc) 


= 1 


4 

1 1 




4 (4A^) 

1 1 
4 (4A^) 


ijkl 


ikr 


32 


kkr 


{N^ - 1 )^ 


(B2) 


where we have used (lAlOb and (lAl lb to get to the second line, 
and the fact that = 0 in the penultimate line, which is a 
consequence of (IA7b and the antisymmetry of f. Using the 
same set of identities, we find that 


T2 — Tr(AiAj)Tr(Ai;A;)Tr(Aj AiAsAr)Tr(ArAsA;Afc) 




fj^S^jSsr 


Ir r 


^p^^sr^kl “t“ gR-rszR'/fcz 


1 

m 


SijSkiSs 

ijklsr 


1 / 1 


4 V 4Af 


= 7 7117 -1) • 


(B3) 


with the coefficients: 


Similarly, one can deduce that 


ro,t 

Uij 


2 If 

- + x Gr,fcG,,fc, 
Us 2 Jk 

f Gp fcGi^fcGj^fe. 

Jk ’ ' ' 


1 

m 


(AI 8 ) 


Trr 



^^m+1 


fN\ 


(B4) 
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